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Abstract: We examine a double trace deformation of SU{N) Yang-Mills theory which, for 
large iV and large volume, is equivalent to unmodified Yang-Mills theory up to 0(1 /N 2 ) 
corrections. In contrast to the unmodified theory, large N volume independence is valid 
in the deformed theory down to arbitrarily small volumes. The double trace deformation 
prevents the spontaneous breaking of center symmetry which would otherwise disrupt large 
./V volume independence in small volumes. For small values of N, if the theory is formulated 
on M 3 x S 1 with a sufficiently small compactification size L, then an analytic treatment of the 
non-perturbative dynamics of the deformed theory is possible. In this regime, we show that 
the deformed Yang-Mills theory has a mass gap and exhibits linear confinement. Increasing 
the circumference L or number of colors N decreases the separation of scales on which the 
analytic treatment relies. However, there are no order parameters which distinguish the small 
and large radius regimes. Consequently, for small iV the deformed theory provides a novel 
example of a locally four-dimensional pure gauge theory in which one has analytic control 
over confinement, while for large iV it provides a simple fully reduced model for Yang-Mills 
theory. The construction is easily generalized to QCD and other QCD-like theories. 
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1. Large iV volume independence 

The large N limit of SU (N) Yang-Mills theories, when formulated on toroidal compactifi- 
cations of M d , are independent of volume provided the (Zjv) rf center symmetry is not spon- 
taneously broken [1, 2, 3]. 1 However, above two dimensions, center symmetry does break 
spontaneously when the (smallest) compactification circumference L is less than a critical 
size L c [4]. (If just one dimension is compactified, then this center-symmetry breaking transi- 
tion is the usual thermally induced deconfinement transition.) In four dimensions, the critical 
size L c is approximately A -1 where A is the MS strong scale of the theory [4, 5]. 

Notwithstanding the limitation to L > L c , the volume independence of large N Yang- 
Mills theory ("partial reduction") has practical utility for lattice studies, because simulations 
on lattices of size (L c ) d are sufficient to extract infinite volume properties of large N Yang- 
Mills theory [4, 6, 7, 8] . But it would be even more helpful to have a formulation of the theory 
in which volume independence holds for arbitrarily small volumes — since this allows one to 
reduce the lattice all the way down to a single site. 

1 Center symmetry transformations are gauge transformations which are periodic only up to an element of 
the center of the gauge group. Volume independence applies to the leading large N behavior of expectation 
values and connected correlators of topologically trivial Wilson loops. 
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Several schemes for preserving volume independence in arbitrarily small volumes have 
been proposed. In so-called quenched reduced models, one constrains the eigenvalues of link 
variables (or Wilson lines) in a manner which prevents Wilson lines from acquiring expectation 
values [3]. 2 In the N — > oo limit, quenched reduced models correctly reproduce properties of 
infinite volume Yang- Mills theory. However, corrections to the N = oo limit scale as 1/N, not 
1/N 2 , in quenched reduced models. This makes extracting large N properties from numerical 
simulations of quenched reduced models quite challenging. An alternative proposal, known 
as twisted reduced models, involves modifying the Wilson action of a single-site model so 
that the action explicitly disfavors configurations in which Wilson lines in different directions 
mutually commute [10, 11]. Unfortunately, this clever scheme fails to work sufficiently close 
to the continuum limit [12, 13]. In essence, the penalty imposed by the twisting of the action 
is insufficient to overcome entropic effects which favor breaking of the center symmetry. 

If light adjoint representation fermions are added to an SU(N) Yang-Mills theory, and 
periodic (not anti-periodic) boundary conditions imposed on the fermions, then the fermion 
contribution to the Wilson line effective potential stabilizes the unbroken center symmetry 
phase. Hence these QCD-like theories satisfy large-TV volume independence for arbitrarily 
small volumes [14]. (In addition, in the large- N limit C-even observables coincide between 
these theories and corresponding theories, in sufficiently large volume, with fermions in the 
rank-two symmetric or antisymmetric tensor representations [15, 16, 17].) 

Motivated by this fermion-induced stabilization of center symmetry, in this paper we 
introduce a simple scheme for preserving volume independence in pure Yang-Mills theory. 
We add double trace terms to the action which prevent spontaneous breaking of center sym- 
metry, while simultaneously perturbing the dynamics of the unbroken symmetry phase only 
by 0(1/ N 2 ) corrections. This leads to a "stabilized reduced model" which reproduces the 
dynamics of infinite volume Yang-Mills theory up to corrections which scale as 1/N 2 . The 
construction may be easily generalized to other QCD-like theories with matter fields in rank- 
one or rank-two representations. 

In addition to providing a simple large N reduced model, the deformed Yang-Mills theory 
is interesting in its own right when N is not large. 3 When formulated on M 3 x S , we show 
that the large distance dynamics of the theory is analytically tractable provided NAL <C 1 • 
In this regime, a semiclassical analysis (closely related to Polyakov's classic treatment of 3d 
SU(2) adjoint Higgs theory [22]) reveals the existence of a mass gap and area law behavior of 
spatial Wilson loops. It is noteworthy that our compactified, deformed Yang-Mills theory is 
an analytically tractable confining theory with no fundamental scalar fields or supersymmetry, 
in contrast to other instructive models of confinement [22, 25, 26, 27]. The confinement mech- 
anism involves the formation of a dilute plasma of magnetic monopoles (and antimonopoles) 
carrying topological charge ±l/iV. 4 

2 See Ref. [9] for an extended discussion of quenched and twisted reduced models. 
3 See also related recent work in Refs. [18, 19, 20, 21]. 

Confinement due to such topological objects has been previously discussed in, for example, Refs. [23, 24] 
and references therein. What is novel about our deformed theory at small iVAL is that this confinement 
mechanism operates in a regime in which one has analytic control over the long distance dynamics. 
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2. Deformed Yang-Mills theory 



We consider pure Yang-Mills (YM) theory with gauge group SU(N) defined on the four 
manifold M 3 x S 1 , with the S 1 having circumference L. The extension to multiple compactified 
dimensions will be discussed below, but we begin with a single compactified dimension to 
simplify the exposition. We start with the usual continuum action, 



or a lattice formulation with the Wilson action, 

S YM = ^ £tr(E/[0p] + E/[0p]t), (2.2) 

peA 4 

where the sum is over all oriented plaquettes, A4 is the four dimensional spacetime lattice, and 
U[dp] denotes the usual product of link matrices around the boundary of plaquette p. The 
lattice coupling (3 = 2/g 2 . In our discussion, we will use both continuum and lattice formu- 
lations, and benefit from both perspectives. As usual, a key virtue of the lattice formulation 
is that it provides an explicit non-perturbative definition of the theory. 

Let fi(x) = V ( e l S dx i A i(*> x i)^ denote the Wilson line (or Polyakov loop) operator — the 
holonomy of the gauge field around a circle wrapping the S 1 and sitting at the point x G R 3 . 
We will construct a deformation of the Yang-Mills action on our compactified geometry by 
adding terms, respecting all symmetries of the unmodified theory, built from the Wilson line 
operator. The deformed action is given by 

^deformed = gYM + A g ? ^3) 

with 

AS^ [ ^P[fi(x)] (2.4a) 

in the continuum, or 

AS = a73 E P I n ( x )l ( 2 - 4b ) 

* xGA 3 

on the lattice. In the lattice form, N t = L/a denotes the size of the lattice in the compact- 
ified direction and A3 C A4 is a three dimensional sublattice of the four dimensional lattice 
corresponding to a fixed Euclidean time-slice. We want the deformation potential to 
guarantee the stability of the phase with unbroken center symmetry (at small volume). It 
will be chosen to have the form 

[N/2] 

P[n}= a n \tr(n n )\ 2 , (2.5) 

with positive coefficients {a n } (and [N/2\ denoting the integer part of N/2). In other words, 

is a sum of the double trace operators tr(f2 n ) tv(fl n y . When considering the large 
limit, the coefficients {a n } will be held fixed as Af — ► 00. 
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Under a center symmetry transformation by some element z £ the Wilson loop 
tr(f2 p ) is multiplied by z p . The value of p mod N, which determines the Z^r representation, 
is referred to as the iV-ality. It will be important that the deformation potential (2.5) only 
contains, by construction, absolute squares of Wilson loops with non-zero iV-ality. 

If P[fi] were simply proportional to |tr 0| 2 , with a sufficiently large positive coefficient, 
then this would prevent breaking of the center symmetry with (tr(fi)) as an order parameter. 
But if N > 3 then this single term is not sufficient to prevent any spontaneous breaking 
of center symmetry, as this term alone does nothing to prevent tr(Q 2 ) from developing an 
expectation value. In other words, a stabilizing term proportional to |tr(S7)| 2 cannot prevent 
Z^r breaking to Z2 (assuming N is even) with tr(0 2 ) as an order parameter. Adding an 
additional stabilizing term proportional to |tr(f2 2 )| 2 could prevent such a breaking to Z2, but 
does not prevent breaking to Z3 (if N mod 3 = 0), or to any larger discrete subgroup of Z^v- 
This is why we have allowed P[£i\ to include terms up to |tr 0^/2] p_ 

We will argue that the deformed theory satisfies the following: 

i) For suitable choices of the deformation parameters a n (i.e., each coefficient sufficiently 
large and positive) the stabilizing potential (2.5) will prevent the Zjy center symmetry 
from breaking to any subgroup. 

ii) In the N — > 00 limit, pure Yang-Mills theory on IR 4 is equivalent to the deformed theory 
formulated on any R 3 x S 1 (for choices of the {a n } satisfying point i). This equivalence 
applies to expectation values of Wilson loops (on IR 4 ), or the leading large N behavior 
of their connected correlators. In the lattice formulation, the number of sites in the 
compactified direction may be reduced to one. 

Hi) When NAL -C 1 limit, the deformed Yang-Mills theory is solvable in the same sense 
as the Polyakov model. The existence of a mass gap and linear confinement can be 
shown analytically. One can regard this regime as having spontaneous breaking of the 
SU(N) gauge symmetry down to [/(l)^" 1 , but this is a perturbative gauge dependent 
description with no well-defined invariant content. 

iv) There exist no order parameters which can distinguish the NAL <C 1 "Higgs" regime 
from the NAL 1 regime in which gauge symmetry is "restored" . 

As noted earlier, pure Yang-Mills theory on IR 3 x S 1 satisfies volume independence (in 
the large N limit) so long as the Z^r center symmetry remains unbroken [2, 3, 4]. We will 
show that this is equally true for the deformed theory; the additional terms in the action 
do not affect the proof that volume independence, at N = 00, is an automatic consequence 
of unbroken center (and translation) symmetry. In the undeformed theory, the unbroken 
center symmetry phase is the low temperature confined phase, L > L c with L c ~ A -1 the 
inverse deconfinement temperature. But in the deformed Yang-Mills theory, the stability of 
the unbroken center symmetry phase is enforced by hand for all values of L. 

Large N volume independence, and the large iV equivalence between ordinary Yang- 
Mills theory on R 4 and the deformed theory on IR 3 x S 1 , may be demonstrated by comparing 
Dyson-Schwinger equations (i.e., Migdal-Makeenko loop equations) for expectation values 
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ordinary Yang-Mills 



deformed Yang-Mills 




Figure 1: Large N equivalences relating ordinary and deformed SU(N) Yang-Mills theories, as a 
function of the size L of the periodic volume in which the theories are defined. In deformed YM 
volume independence holds for all L, while in ordinary YM volume independence fails below a critical 
size, L < L c , (shaded region) due to spontaneous breaking of center symmetry. This prevents reduction 
all the way down to a single-site matrix model for ordinary YM. Large N equivalence holds between 
ordinary and deformed YM theories as long as center symmetry is unbroken in ordinary YM. Large N 
volume independence is a type of orbifold equivalence, as discussed in Ref. [14]. The combination of 
volume changing orbifold projections in the deformed theory, along with the deformation equivalence 
in sufficiently large volume, provides a useful equivalence between deformed YM in small volume 
and ordinary YM in large volumes. In particular, a single-site matrix model of the deformed theory 
will reproduce properties of ordinary Yang-Mills theory in infinite volume. The construction can be 
generalized to QCD, with the deformed theory providing a fully reduced matrix model for QCD. 



and correlators of Wilson loops, or alternatively by comparing the large N classical dynamics 
that may be derived by using appropriate large TV" coherent states [2, 30, 31, 29]. Figure 1 
summarizes the relation between the large N limits of ordinary and deformed Yang-Mills 
theories. As long as the Z^r center symmetry is not spontaneously broken, the dynamics of 
the theories denned by S YM and 5' deformed a re indistinguishable at leading order in the 1/N 
expansion. In particular, glueball spectra of the two theories can differ only by order 1 /N 2 
effects. Agreement up to 0(1/N 2 ) terms also applies to the string tension which characterizes 
the area law behavior of large Wilson loops. 

The large N equivalence between ordinary Yang-Mills theory and the deformed theory 
in large volume, combined with the large N volume independence of the deformed theory, 
circumvents the problems with previous formulations of reduced models for pure Yang-Mills 
theory. Unlike the original Eguchi-Kawai model [1], its twisted variant [10, 11], and the partial 
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reduction of Refs. [4, 5], the equivalence to deformed YM theory remains valid in the limit 
of zero compactification radius, irrespective of the value of the (bare) gauge coupling. And 
for finite N, corrections to the large N limit scale as 1/N 2 , not 1/jV as in quenched reduced 
models [3]. Consequently, it should be possible to study the deformed theory, for relatively 
modest values of iV and vanishingly small volume, and obtain accurate results for properties 
of ordinary Yang-Mills theory, in the large N limit, on M 4 . As we discuss below, it is also 
instructive to study the deformed theory for small values of N. In this regime, it will be seen 
to provide a novel example of a confining theory, only involving an SU(N) gauge field, which 
is analytically soluble. 

2.1 Stabilization of center symmetry 

The possibility of preventing spontaneous breaking of center symmetry through the addition 
of a deformation potential of the form (2.5) is largely self-evident. A positive coefficient a n 
suppresses configurations in which tr(J7 n ) is non-zero. Although the pure-gauge dynamics 
of the undeformed theory, in small volume, leads to an effective potential for the Wilson 
line which is minimized when is an element of Z^r, adding the deformation potential P[£l] 
changes the shape of the Wilson line effective potential. For sufficiently large values of the co- 
efficients {a n }, the effective potential will be minimized by configurations in which tr(O n ) = 
for 1 < n < \_N/2\ . This implies that tr(J7 n ) = for any integer n which is non-zero mod- 
ulo N because, for 5i7(iV)-valued matrices, tr(J7 n ) is not independent of lower order traces 
when n > \_N / 2 J . Vanishing of these traces implies that the eigenvalues of 0, are uniformly 
spaced around the unit circle, so that the set of eigenvalues is invariant under Zjv transfor- 
mations (which multiply every eigenvalue by e 2m l N ). This shows that the center symmetry 
is not spontaneously broken. Henceforth, we assume that the coefficients {a n } of the defor- 
mation potential P[£i\ are suitably chosen so as to enforce unbroken center symmetry for all 
compactification radii. 

This argument may be made much more explicit if one considers small compactifica- 
tions, L -C A , so that (due to asymptotic freedom) the gauge coupling at the scale of the 
compactification is small and the theory is amenable to a perturbative treatment. Quantum 
fluctuations generate a nontrivial potential for the Wilson line [32]. In ordinary Yang-Mills 
theory, integrating out the gauge field (and Faddeev-Popov ghosts) produces the functional 
determinants 



where det+ denotes a determinant in the space of periodic functions with period L. Therefore, 
the effective potential for the Wilson line is 



For constant (or slowly varying) configurations, the evaluation of the functional determinant 
is straightforward and yields [32] 



[det + (-i^, V)]" 1/2 [det + (-^dj)] = [det+C-D^)] 



-i 



(2.6) 



V[n] = L" 1 In det+(-l£y) . 



(2.7) 




(2.8) 



with 
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2 °° 1 

V[n] = --j£-j |tr^-| 2 . (2.9) 
7r z n 4 
n=l 

For sufficiently small L, corrections to this one-loop result are negligible. The effective po- 
tential (2.9) is minimized when the Wilson line is an element of the center, O = e 2mk / N ; 
k = 0, • • • , N— 1, demonstrating the spontaneous breaking of Zjy symmetry, in ordinary 
Yang- Mills theory, for sufficiently small compactifications. 

To force unbroken center symmetry in the deformed theory, the deformation potential 
P[fi] must overcome the effect of the one-loop potential (2.9). A simple specific choice for 
the deformation coefficients which, in the continuum limit and for sufficiently large N, ac- 
complishes this is a n = 4/(-7r 2 n 4 ). For this choice, the deformation potential J* K3 L~ 4 P[0] is 
minus twice the one-loop Wilson line effective potential (2.9) of the undeformed theory, so the 
net effect of the deformation is to flip the sign of the effective potential for the Wilson line. 
The resulting combined potential is minimized when trO n = for all n which are non-zero 
modulo N, indicating unbroken center symmetry. 

2.2 Multiple compactified dimensions 

Instead of compactifying just one dimension, one may consider the gauge theory on R 4 - d xT d , 
where T d is a d > 1 dimensional torus. For simplicity, we will discuss the case of a symmetric 
torus having size L in each dimension. Two particularly interesting cases are d = 3, where 
the L — ► limit will reduce to SU(N) matrix quantum mechanics, and d = 4, where the zero 
size limit will reduce to a matrix model. 

With multiple compactified dimensions, the center symmetry is (Z7v) d . The preceding 
discussion generalizes in a straightforward fashion, but there is one significant change: with 
multiple compactified dimensions there is a richer set of possible symmetry realizations in 
which the center symmetry is partially broken. One needs to craft the deformation potential 
in such a fashion that all unwanted symmetry realizations are suppressed. 

With multiple compactified directions, one may define distinct Wilson line operators 
f2i, 0,2, £ld each of which wraps one elementary cycle of the torus. The classical Yang- 
Mills action suppresses configurations in which Wilson lines in different directions are non- 
commuting. (In other words, configurations with vanishing field strength must have Wilson 
lines which commute and are covariantly constant.) A deformation potential with the form 
(2.5), involving just one Wilson line fij, can prevent symmetry realizations in which the 
Zjv center symmetry associated with direction i spontaneously breaks in a manner which 
is independent of the Z^v symmetries associated with other directions. However, it is also 
possible for the center symmetry to break in a manner which leaves a unbroken subgroup 
involving correlated symmetry transformations in different directions. With, for example, 
d = 2, suppose that the (Z^) 2 symmetry breaks to the diagonal Z^r subgroup. Such a 
symmetry realization does not allow expectations of tr fij or tr f2| to be non-vanishing (for 
any k which is non-zero modulo N). But this realization does allow trfii^ to have a non- 
zero expectation value. To prevent such a symmetry realization, it may be necessary to add 
a |tr SI1O2I 2 term to the deformation potential. 



-7- 



To see that this concern is not moot, it is instructive to examine the one-loop effective 
potential for Wilson lines when there are multiple compactified directions. Once again, if the 
physical size of the ci-torus is much smaller than the inverse strong scale A -1 , then "Kaluza- 
Klein" modes with momenta at and above the compactification scale can be integrated out 
perturbatively. This produces a one loop effective potential for mutually commuting Wilson 
lines given by [33, 34] 

v^,...,^-- yi yy — (2 - io) 

(n 1 ,-,n d )e(Z' 1 -0) \ U l~ l l ~ U d) 



This effective potential has N d degenerate global minima at which each Wilson line is an 
independent element of the center. This demonstrates the spontaneous breaking of {1*N) d 
symmetry (down to the identity) in ordinary Yang-Mills theory on R 4 ~ d x (S l ) d for sufficiently 
small compactifications. 5 Adding a deformation to the Yang-Mills action which is the direct 
generalization of the previously discussed deformation (2.4), namely 



with P[fi] having the form (2.5), will deform the effective potential and can prevent this 
symmetry realization in which the (Zjv) d symmetry breaks down to nothing. 

However, the fact that all terms in the effective potential (2.10) have negative coefficients 
demonstrates that the center symmetry will partially break spontaneously even in the pres- 
ence of stabilizing terms of the form (2.11). In particular, even when all traces wrapping a 
single compactified direction vanish, the potential (2.10) favors configurations in which traces 
wrapping multiple cycles (such as trOiO^) are non-zero relative to the center-symmetry pre- 
serving configuration we are trying to stabilize (for which tr Q™ 1 ■ ■ ■ = for all rn that 
are non-zero modulo N). 

Therefore, to prevent any spontaneous breaking of the center symmetry when there are 
multiple compactified directions one must allow the deformation of the action to include 
absolute squares of order parameters which wrap multiple cycles of the torus. In other words 
we need to have 

AS= f -^P[n u ... ,U d ], (2.12) 

5 More precisely, this demonstrates the spontaneous breaking of (Zn) 4 symmetry down to the identity 
when d = 1 or 2. If d > 2 then there are fewer than two non-compactified directions. In this case, the center 
symmetry cannot break spontaneously for any finite value of N, due to the non-zero probability of fluctuations 
(or tunneling events) which effectively average over all degenerate minima of the Wilson line effective potential. 
However, the probability of such fluctuations vanish exponentially with increasing N. Consequently, in the 
N — > oo limit one can have spontaneous breaking of the (Zjv) d center symmetry even when all directions are 
compactified. In order to have a large N equivalence between the compactified theory and Yang-Mills theory 
on M 4 , one must prevent spontaneous breaking of the center symmetry in the large N limit. Therefore, the 
more complicated form of the stabilizing potential discussed below is necessary for all d > 1. 
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where the deformation potential has the form 

[_N/2\ 

P[Sli,.-. ,Q d ]= Yl a ^-n d |tr (^•••0^)| 2 , (2.13) 

ni,- ,n d =-[N/2] 

with coefficients a ni ... nd which are sufficiently large and positive. 

A simple specific choice for the coefficients which, in the continuum limit and for suf- 
ficiently large N, accomplishes this is a ni ... nd = \ (n\ + • • • + n^) . For this choice, the 
deformation potential is minus twice the one-loop effective potential (2.10) of the undeformed 
theory, so once again the net effect of the deformation is to flip the sign of the effective po- 
tential for the Wilson line. The resulting combined potential has a unique center-symmetry 
preserving minimum. 6 

In the following sections we will, for simplicity, largely focus on the case of just one 
compactified direction. 

2.3 Large N equivalence between ordinary and deformed YM 

The most direct way to demonstrate equivalence between ordinary Yang-Mills theory and our 
deformed theory is to compare the Schwinger-Dyson (or loop) equations for gauge invariant 
observables. As usual, for a rigorous treatment it is appropriate (and convenient) to work with 
lattice regulated formulations of both theories. It is also convenient to consider U{N) gauge 
theories instead of SU(N); the difference in gauge groups only affects subleading 0{1/N 2 ) 
relative corrections to Wilson loop expectation values or connected correlators. 

Let 81 denote an operator which varies individual link fields according to 8g(U[£']) = 
6£e't a U[£], where {t a } is a set of U(N) Lie algebra basis matrices satisfying tr t a t b = \S ab . 
Invariance of the Haar measure implies that the integral of any variation vanishes, 

J dn 51 (anything) = , (2.14) 

where d/z = \\ & dU[£']. Choose (anything) to be e s 6% W[C\, where W[C] = ^tr£/[C] is the 
Wilson loop around some closed contour C. Summing over the Lie algebra index a and the 
link £ yields 

J d^o e s (SS ■ 5W[C] + 5 2 W[C}) = , (2.15) 

6 In the case of R x T 3 , the undeformed theory was studied in detail by Liischer, van Baal and others in 
the mid 80's. (See the reviews [35, 36] and references therein.) The goal of this "QCD in a box" program 
was to use asymptotic freedom combined with the absence of the phase transitions in finite volumes to extract 
lessons about QCD on R 4 . This approach confines the theory to its short distance perturbative regime where 
it is perturbatively solvable. However, features of this small universe are very different from QCD on R 4 . In 
particular, when the physical box size becomes comparable to the inverse strong scale A -1 there is a cross-over 
from a world of hadrons to a world of quarks and gluons. As noted in footnote 5, this cross-over becomes a 
sharp phase transition in the large N limit. 

In our deformed Yang-Mills theory, as will be discussed below, there are no phase transitions as a function 
of box size even in the N = oo limit. For the deformed theory we will find that it is not 1/L which acts 
as an infrared cutoff, but rather 1/LN. If LNA 3> 1 then the deformed theory reproduces the dynamics of 
Yang- Mills theory on M 4 . 
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where the dot product is shorthand for the sum over a and I. Dividing by the partition 
function Z = J cL^lq e s yields relations among expectation values, 

(6S ■ 5W[C}) + (5 2 W[C}) = 0. (2.16) 

These are Schwinger-Dyson equations for Wilson loop expectation values. 

After working out the action of the variations, the result may be expressed in a purely 
geometric form. For lattice gauge theory with the Wilson action one finds [37], 

{W[(dp)C]) - {W[(dp)C]) 

ZCC p\£cdp 

+ ^ t(W[C']W[C"}) . (2.17) 

self— intersections 

Here \C\ is the length of the loop C (i.e., the number of links in the loop), W[(5p)C] denotes a 
Wilson loop which goes around the boundary of plaquette p (which contains a link contained 
in the contour C) and then around the contour C, and dp denotes the oppositely oriented 
plaquette boundary. The sum over self-intersections runs over all ways of decomposing a 
loop C which multiply traverses some link £ into two separate loops, C = C'C" , with the 
associated sign determined by whether C and C" traverse the link t in the same or opposite 
directions. See Ref. [28] for more detailed discussion. 

In the large N limit, with = = f held fixed (where A = g 2 N is the 't Hooft coupling), 
all N dependence disappears. Fluctuations in the values of Wilson loops vanish in this limit 
(their distributions become arbitrarily sharply peaked). This is a reflection of the classical 
nature of the large N limit [2] , and implies that the expectation value of a product of loops 
factorizes, up to 1/N 2 corrections, 

(W[C')W[C"}) = (W[C'))(W[C"]) + 0(1/N 2 ) . (2.18) 

(The 0(1/ N 2 ) remainder is the connected correlator.) Consequently, in the large N limit 
Wilson loop expectation values satisfy a closed set of nonlinear algebraic equations, 

\\C\{W[C)) = Y, E I [{W[(dp)C)) - {W[(dp)C)) 

£cCp\lcdp 

+ E T(W[C']) (W[C"]) . (2.19) 

self —intersections 

The loop equations in the concise form (2.16) are equally valid for the deformed theory. 
The only difference is that S now includes the double trace deformation AS, and this generates 
new terms in the loop equations given by (S(AS) ■ 5W[C]). Just as the usual Wilson action 
leads to terms in which a plaquette is inserted into the loop C, the piece of AS proportional 
to |trO fc | 2 generates terms in the loop equation in which the topologically non-trivial loop 
fi fc (or its inverse) is "sewn" into the loop C (if C contains links pointing in the compactified 



l\C\(W[C}) =E E 
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direction). But because A5 contains absolute squares of traces, each such term is multiplied 
by the complex conjugate of the trace of the inserted loop. Hence, 

(s(as) ■ 5W[c\) = h[c\ (w[n k c\ w[n- k \) , (2.20) 

where Q fc C denotes a loop obtained by concatenating fi fc and C at their intersection links, 
and the coefficients b^[C] are proportional to a\u but also depend on the number of links in 
C which point in the compactified direction. The essential point is that the variation acts 
on one of the two traces comprising the double trace deformation, leaving the other trace 
unchanged. In the large N limit, due to factorization, 

(w[n k c}w[n- k \) = (w[n k c\) (w[n- k \) + o(i/n 2 ) . (2.21) 

But for any k which is non-zero modulo N, W[Q~ k ] transforms non-trivially (acquiring a 
phase e -' 2mk / N 'j under a Z^r center symmetry transformation. Hence its expectation value is 
an order parameter for the center symmetry and (VF[f2 -fc ]) must vanish in any phase with 
unbroken center symmetry. As discussed above, the deformed theory is constructed so as to 
ensure unbroken center symmetry for all compactification radii. Consequently, all additional 
terms in the loop equations generated by the deformation of the action vanish in the large TV" 
limit, 

(6(AS) ■ 5W[C] ) = 0(l/N 2 ), (2.22) 

implying that Wilson loops in the original and deformed Yang-Mills theory satisfy identical 
large-iV Schwinger-Dyson equations. 

Ordinary Yang-Mills theory has unbroken center symmetry only for sufficiently large 
compactifications, L > L c . The coinciding large N loop equations in ordinary and deformed 
Yang-Mills theories imply that Wilson loop expectation values in these two theories have 
identical large N limits when L > L c . 7 

The same approach may be used to compare the Schwinger-Dyson equations satisfied by 
connected correlators of two or more Wilson loops, with exactly the same conclusion: the 
leading large N behavior of connected correlators coincide between ordinary and deformed 
Yang- Mills theories, provided L > L c . Thus, in sufficiently large volume the net effect of 
the double trace deformation on the dynamics of the theory is 0(1/N 2 ), and vanishes in the 
large iV limit. In another words, the physics of the deformed Yang-Mills theory depends on 
the deformation parameters {a^} only in the combination ai/N 2 which vanishes at N = oo. 
This demonstrates the nonperturbative equivalence of ordinary Yang-Mills theory and the 
deformed YM theory, formulated on M 3 x S 1 (or more generally, on any toroidal compacti- 
fication of flat space), provided the compactification size is above the critical size for center 
symmetry breaking in the undeformed theory. 

7 This argument, that coinciding loop equations imply coinciding expectation values, glosses over the possi- 
bility that the infinite set of loop equations may have multiple solutions which respect center symmetry, with 
different theories potentially corresponding to different solutions of the same set of equations. The alternative 
approach of comparing the N — oo classical dynamics generated by appropriate large N coherent states, dis- 
cussed in Ref. [29], eliminates this loophole and demonstrates equivalence in any phase of the theories which 
satisfy the necessary and sufficient symmetry realization conditions. 
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2.4 Large iV volume independence of deformed YM theory 

Unbroken center symmetry is necessary and sufficient for the validity of the large N volume 
independence of Yang- Mills theory (or more general gauge theories containing adjoint repre- 
sentation matter fields). This may be demonstrated by comparing large N loop equations, 
or the N = oo classical dynamics generated by suitable coherent states [14]. Corrections to 
this equivalence for finite iV scale as 1/N 2 . The loop equation analysis is very similar to that 
sketched above. In the large N loop equations for topologically trivial Wilson loops, one finds 
that the only volume-dependent terms (arising from self-intersections) automatically vanish 
as long as the center symmetry is not spontaneously broken. 

The analysis of large N volume independence in Ref. [14] applies equally well to the 
deformed theory which, by construction, has unbroken center symmetry for any compactifi- 
cation radius. Because the double trace operators in P[Cl] are squares of loops with non-zero 
Y-ality, the presence of the deformation potential P[fi] has no effect on the large N classical 
dynamics of center-symmetry symmetric states. Consequently, deformed Yang-Mills theory, 
in the large N limit, is completely volume independent. 

In the lattice formulation, if one compactifies all directions then one may reduce the 
lattice size all the way down to a single site, in which case the theory becomes a simple 
matrix model of Wilson lines running in each lattice direction with action, 

5s1n"e = f E tr (WjSiln] + o,o,(M):) + p[fi ls . . . , n d] . (2.23) 

i>j=i 

The large N limit of this matrix model will reproduce the leading large N behavior of expec- 
tation values and connected correlators of Wilson loops in uncompactified Yang-Mills theory. 
As discussed in the Introduction, the single-site deformed Yang-Mills theory (2.23) provides 
a simple generalization of Eguchi-Kawai reduction which is valid for any value of the lattice 
coupling p. 

2.5 Addition of matter fields 

Consider adding TVf species of matter fields (either fermions or scalars) in the fundamental 
representation to SU(N) Yang-Mills theory, either ordinary or deformed, with one dimension 
compactified. The addition of fundamental representation matter explicitly breaks the Zjy 
center symmetry. However, if TVf is held fixed as N — > oo, then the fundamental representation 
matter fields have only a subleading 0(Nf/N) effect on the gauge field dynamics. As a result, 
everything discussed above remains valid. That is, the leading large N behavior of expectation 
values or connected correlators of Wilson loops in the undeformed theory, in sufficiently 
large volume, coincide with the corresponding observables in the theory, in arbitrary volume, 
deformed by the addition of the stabilizing potential P[fi]. In addition, one may also show 
that the same equivalence applies to the leading large N behavior of mesonic expectation 
values and connected correlators. 8 Note, however, that these large N equivalences cease to 
apply if Nf/N is held fixed as N — > oo [14]. 

8 This is easiest to understand by considering the equivalent gluonic observables produced by integrating 
out the matter fields. For the case of large N volume independence in the undeformed theory, see Ref. [14] for 
details. The presence of the deformation potential does not affect this analysis. 
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If adjoint representation matter fields are added to the theory (ordinary or deformed), 
then the large N equivalences discussed above also remain valid. Adding adjoint representa- 
tion fields enlarges the natural set of gauge invariant observables from simple Wilson loops to 
Wilson loops decorated by arbitrary numbers of insertions of adjoint matter fields. But the 
presence of adjoint matter fields preserves the center symmetry of the underlying Yang-Mills 
theory. As a result, the above-described comparison of large N loop equations (or large N 
classical dynamics) between the ordinary and deformed theories immediately generalizes to 
the case of Yang-Mills theories with adjoint matter, with exactly the same conclusions. 9 

Finally, one may also consider the addition of matter fields in rank-two antisymmetric or 
symmetric tensor representations (yielding theories we will refer to as QCD(AS) or QCD(S), 
respectively). The presence of fields in these representations reduces the U(l) center sym- 
metry of U(N) Yang-Mills theory down to Z2. Given the central role the center symmetry 
played in the above large N equivalences, one might think this reduction in center symmetry 
would destroy these large TV" equivalences. This is not the case. One way to see this is to 
note, as discussed in Ref. [14], that volume-dependent terms in the N = 00 loop equations 
only appear if loops with non-zero winding number around the compactified direction acquire 
non-zero expectation values. The addition of the deformation potential P[fi] prevents topo- 
logically non-trivial Wilson loops from acquiring non-zero large N expectation values, even in 
small volumes. Consequently, the situation is analogous to the pure gauge case: the leading 
large N behavior of expectation values or connected correlators of single trace observables in 
QCD(AS/S) in sufficiently large volume coincide with the corresponding observables in the 
theory modified by the addition of the deformation potential in arbitrary volume, as depicted 
in Fig. 1. 10 

3. Confinement at small radius and small N 

When compactified on a small circle, L <C A -1 , the gauge coupling of the deformed theory is 
small at the compactification scale, g 2 (l/L) <C 1. As discussed earlier, the combined potential 

9 For a detailed discussion of loop equations in theories with adjoint matter, see Ref. [28]. 

10 Another way to understand this is to note the existence of a large N equivalence (so-called "orientifold 
equivalence") between theories with rank-two symmetric or antisymmetric representation matter and corre- 
sponding theories with adjoint representation matter ["QCD(adj)"] [15, 16, 17]. This large N equivalence 
applies to the charge-conjugation even sectors of the two theories, and only holds if charge conjugation sym- 
metry is not spontaneously broken. When, for example, the matter fields are fermions with periodic boundary 
conditions, examination of the Wilson line effective potential shows that QCD(AS/S) does spontaneously 
break both charge conjugation and center symmetry when compactified with sufficiently small size [15]. But 
the addition of a deformation potential of the form (2.5) (with sufficiently positive coefficients) will prevent this 
spontaneous symmetry breaking, just as it does in the pure Yang-Mills case. Since QCD(adj) satisfies large N 
volume independence (as long as its center symmetry is not spontaneously broken), the same large N volume 
independence must also apply to QCD(AS/S) (as long as charge conjugation is not broken). In sufficiently large 
volumes, there is no reason to believe that charge conjugation symmetry breaks spontaneously in QCD-like 
theories with rank-two tensor representation matter. Therefore, large N orientifold equivalence combines with 
large N volume independence of QCD(adj) to imply volume independence in QCD(AS/S) as long as center 
and charge conjugation symmetries are not spontaneously broken — which is what the deformation potential 
ensures. 
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V[0] + P[£l] is minimized when 



U = Diag ( 1, e 2 ^ N , e 4 ™^, ■■■ , e 2 ^^^ ) , (3.1 



up to conjugation by an arbitrary SU(N) matrix. Working in a gauge in which 0, is diagonal, 
and using (temporarily) gauge-dependent language, this configuration may be regarded as 
breaking the gauge symmetry down to the maximal Abelian subgroup, 

SU(N) U(1)( N -V . (3.2) 

Modes of the diagonal components of the SU (N) gauge field with no momentum along the 
compactified £4 direction describe photons associated with the Cartan subgroup of SU(N). 
Modes of the diagonal components of the gauge field with non-zero momentum in the com- 
pactified direction form a Kaluza-Klein tower and receive masses which are integer multiples 
of 2ir/L. The off-diagonal components of the SU(N) gauge field describe Kaluza-Klein tow- 
ers of W^-bosons which are charged under the unbroken [/(l)^ -1 gauge group. The non-zero 
value of A4 = —{i/L) In $7 shifts the masses of these off-diagonal components by multiples of 
2n/(NL). The net effect is that there are charged VF-bosons with masses 

mWk = i2 WL J k=l,2,---,. (3.3) 
For later convenience, we define mw to be the mass of the lightest W bosons, 

mw = m . (3.4) 

This is the mass scale below which the dynamics is effectively Abelian. 11 Note that, at fixed 
L, the lightest W bosons have masses which become small when N — ► 00. This will be 
important in the discussion of the large N behavior of the deformed YM theory. But first, in 
this section, we consider the dynamics of the deformed theory when N is fixed and small. 

The N— 1 photons of the Cartan subgroup do not couple (directly) to the Wilson line and 
remain massless to all orders in perturbation theory. Thus, a strictly perturbative analysis 
would lead one to expect that the deformed theory, for sufficiently small L, would have a 
non-confining Coulomb phase. We will see that this is incorrect — nonperturbative effects 
lead to the generation of a mass gap and produce confining long distance physics. 

The analysis of non-perturbative properties in our compactified deformed Yang-Mills 
theory is very similar to Polyakov's treatment of the 3d Georgi-Glashow model [22]. But 
instead of a three-dimensional theory with a non-compact Higgs field, we have a compactified 
four-dimensional theory with the group-valued Wilson line O serving as a compact Higgs 
field. For theories involving massless complex fermions, the difference between compact and 



11 Fluctuations in the eigenvalues of 51 away from the minimum (3.1) correspond to neutral "Higgs bosons". 
The masses of these fluctuations depend on the coefficients {a n } of the deformation potential but parametrically 
arc of order y/\/L ~ ^XN mw- This scale will be large compared to the mass scale of the non-perturbative 
3d dynamics, and these fluctuations will play no role in the following discussion. 
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noncompact Higgs systems can be major [38]. However, in our case, the differences relative 
to Polyakov's classic discussion are rather minimal. 

Due to the SU(N) — > [/(l)^ -1 gauge symmetry "breaking", there exist topologically 
stable, semiclassical field configurations, namely monopoles [32]. At the center of a monopole, 
one pair of eigenvalues of the Wilson line become degenerate. For fundamental (i.e., minimal 
action) monopoles, this will be a pair of eigenvalues which are nearest-neighbors at infinity. 
If the adjoint Higgs field was noncompact, then there would be N— 1 species of fundamental 
monopoles. This follows from the topological considerations: the second homotopy group 
iT2[SU(N) /U(1) N ~ 1 ] = tti[U(1) n ~ 1 ] = Z^" 1 , implying that fundamental monopoles come in 
N— 1 varieties. However, with a compact Higgs field there is an extra fundamental ("Kaluza- 
Klein") monopole which arises due to the fact that the underlying theory is formulated on a 
cylinder, R 3 x 5 1 , or equivalently that the configuration space of O is compact. 

The monopoles may be characterized by their magnetic charges, topological charge, and 
their action. The magnetic charges of the N different types of fundamental monopoles are 
proportional to the simple roots and affine root of the Lie algebra of the unbroken U(1) N 
gauge group. 12 The simple roots are given by 13 

a\ = (1, —1, 0, . . . , 0) = ei — e.% , (3.5a) 
a-2 = (0,1,-1, ,0) = e 2 - e 3 , (3.5b) 

a N -i = (0, . . . , 0, 1, -1) = e N -i - e N , (3.5c) 

and the affine root is 

N-l 

&N = - a i = ( _1 > 0, 0, . . . , 1) = e/v - ei . (3.6) 
i=i 

For later convenience, let A^ ff denote the affine (extended) root system of the the associated 
Lie algebra, 

A aff - {"i> a 2, • • • ,ajv-i,Q!Jv} • (3.7) 

It is the affine root system which is relevant for compact Yang-Mills Higgs systems. The roots 
ai G A° ff obey 

ai ■ aij = 25ij - Sij+i - Sij-i , i,j = l,...N. (3.8) 

The form (3.8) of these inner products will translate into self and nearest neighbor interactions 
between monopoles in the Dynkin space. The above choice of basis is natural due its visual 

12 Even though the gauge symmetry "breaking" is SU (N) — > [/(l)^ -1 , for ease of presentation it is convenient 
to add an extra photon to the original theory and discuss U(N) — > U(1) N . This simplifies the discussion of 
charge assignments of monopoles, and the affine roots of the associated Lie algebra. In the continuum limit 
of the theory this extra photon completely decouples from the other degrees of freedom and may simply be 
ignored. It should not be confused with the AT— 1 photons which have non-trivial nonperturbative dynamics. 

13 This set of simple roots corresponds to choosing Lie algebra generators normalized to satisfy tr t a t b — 5 ab , 
instead of | 8 ab as in the previous section. 
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simplicity, but the inner products (3.8) of the roots of the associated Lie algebra are basis 
independent. 

Let = ^e^ uM F u \ denote the U(1) N valued 3d magnetic field, with conventional 
perturbative normalization. (In a gauge where O is diagonal, is just the list of diagonal 
elements of the original non-Abelian field strength, multiplied by 1/g.) The magnetic charges 
of a monopole of type i = 1, ■ ■ • , N are given by the root ai (up to a factor of 2n/ 'g), 

/ d~E • F = — on [type (i) monopole] . (3.9) 

Js? 9 

(The S 2 is an arbitrarily large sphere in K 3 . The flux is independent of the value of x A at 
which the integral is performed, as the long distance monopole fields are independent of x 4 .) 

The topological charge is correlated with the magnetic charge of the monopole. For 
fundamental monopoles with magnetic charges ctj G A^ ff , the topological charge is 

J^trF^ = l. (3.10) 
: X 5i lbir z TV 

For antimonopoles with magnetic charges — Qj, the topological charge v = —1/N. 

The electric charges of TV-bosons may also be simply expressed in terms of the affine 
roots A° ff . The lightest W bosons, with mass m\y, may be labeled by a single root which 
gives their electric charges (up to a factor of g), 

Qw ai = ga-i. (3.11) 

TV-bosons in the next heavier multiplet are labeled by a pair of neighboring roots, and have 
charges 

Qw ai+aw = 9 («i + a i+1 ) , (3.12) 

etc. Dot products of the TV-boson charges and monopole charges obey the Dirac quantization 
condition, 



2tt 

Qw a . ■ Qhi a . = goti- — ctj = 2n(25ij - Si 4+i - Sij-i) 

J g 



4tt , for i = j ; 
-2vr, fori=j±l; (3.13) 
, otherwise. 



Conjugation by a 7hfq "shift" matrix, which is part of the global gauge symmetry, cyclically 
permutes the Wilson line eigenvalues and hence cyclically permutes the N different species of 
fundamental monopoles. The presence of this symmetry (which is one of the features which 
distinguishes compact and non-compact Higgs systems) guarantees that the N different types 
of fundamental monopoles have identical values of the action. Monopole solutions are self- 
dual, 

Ffxu = Ffiv , (3.14) 
and hence the Yang-Mills action of a fundamental monopole (or antimonopole) is 



— tr F F^ 

'xS 1 2 9 



8tt 2 , 8tt 2 



9 2 11 g 2 N' 



(3.15) 
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After adding the contributions of the deformation potential and the induced one-loop 
effective potential V[fi], the complete monopole action will differ from this value. But the 
deviation is perturbative in g 2 , so the monopole action 



So = S YM + AS + S 1 _ loop = %^ + 0(1) • 



(3.16) 



The correct infrared description of the deformed Yang-Mills theory on M 3 x S 1 at small 
radius is generated by a dilute gas of monopoles (and antimonopoles) of N different types, 
interacting via the species-dependent long range Coulomb potential, 



/ 2tt\ 2 25ij - Sij+i — Sij-i 



47r|r| 



(3.17) 



(The overall sign is plus for monopole-monopole, and minus for monopole- antimonopole.) 
Hence, we are dealing with a multi-component classical plasma, with nearest-neighbor inter- 
actions in the Dynkin space. As with any classical plasma, this system will exhibit Debye 
screening. The field due to a static external magnetic charge will fall exponentially with 
distance, |F| ~ e~ mDT /r, with m^, 1 the characteristic Debye screening length. This implies 
that external fields cannot propagate coherently over distances large compared to the Debye 
length, which will be the longest correlation length in the system. The Debye mass mo will 
appear as a dynamically generated photon mass. This will be shown explicitly. 

For momentum scales small compared to the lightest W mass, the equilibrium dynamics 
is correctly represented by a grand canonical ensemble of all types of monopoles and anti- 
monopoles. Consider a configuration in which there are monopoles and antimonopoles 



of types i = 1, ■ ■ • ,N, located at positions rjj. , k = 1, ■ ■ ■ , nW and r| , I = 1, ■ 
spectively. The magnetic field generated by this ensemble of magnetic charges is 



, n 



(0 



re- 



N 



B ( x ) = Ev a * 



E 

k=l 



x — r 



(0 



47r|x 



.«|3 



E 



x — r 



(0 



1 47r|x 



,W|3 



(3.18) 



The action of such a monopole configuration is the sum of the monopole self-energies plus 
their potential energy due to Coulomb interactions, 



'monopole— gas 



N 



(3.19) 



with 

-Suit = 

and 



N 



E ai " a i 



„(») n U) n(i) nU) n (i) nU) 



fc=l 1=1 



k=l 1=1 



E E G(t®-tP) + E E G(rW-rW) - 2 £ £ G(rW-f p>) 

fc=i i=i 

(3.20) 

1 



G(r) = 



47r|r| 



(3.21) 
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The grand canonical partition function of this multi-component Coulomb gas is 



^niEsEL ,n *s? /.., n *?> > «-*• . <w 




nW=0 " K fc=l " /K i=l 



where 



C = c e" 50 = A ( 5 2 iV)~ 2 e- A5 e -8*WW) (3 . 23 ) 



is the monopole fugacity. The prefactor C represents the one-loop functional determinant in 
the monopole background. Extracting the zero- modes of the small fluctuation operator via 
the usual collective coordinate procedure leads to factors of (g 2 N)~ 2 niyy . (See the appendix 
for details.) If the coupling is evaluated at the scale mw, which is natural for this problem, 
then the non-zero mode part of the one-loop determinant merely gives rise to an overall 
dimensionless (and N independent) coefficient A. In the final form of (3.23), AS* denotes the 
deformation term in the action (2.4) evaluated in the background of a fundamental monopole. 
This is an 0(1) number, independent of the coupling g 2 , whose explicit value depends, of 
course, on the deformation parameters {a n }. [For large N, AS scales as 0(1/N).] 

Using the fact that G(r) is the Green's function for the 3d Laplacian, this partition 
function can be exactly transformed into a 3d scalar field theory with an iV-component real 
scalar field, 

N 

Z = / Yl Vai e" Sdual[<T] , (3.24) 



N 



where 

5 dual = f r 1 (jL\ 2 (V<T) 2 _ c £ CQs{ai . a) 

J R 3 L2L \2ir/ ^ J 

To verify this, it is easiest to start with the functional integral (3.24), rewrite the cosines in 
terms of exponentials of <r, expand the exponential of each of the resulting interaction terms 
in a power-series in e~ So e ±JQi ' <T , and then perform the functional integral over a. The scalar 
fields <7j appearing in this representation are dual fields for the 3d Abelian gauge fields A 1 ^. 14 
The fields {<7j} should be regarded as compact scalar fields defined modulo 2tt. In addi- 
tion to invariance under 2ir shifts in any component of <x, note that the monopole induced 



14 In three dimensions, Abelian duality relates a photon to a compact scalar. With cr-' (x) the compact scalar 
dual to the photon ^4^(x) of the j'th U(l) subgroup, the Abelian duality relations are 

*da j = \L Im(r) F« , r^- 1 ) = ^ + ± , F$ = e^ p d p a> . 

The 3d Maxwell action becomes j^j (F^) 2 = LI ™ ( J ] (F$) 2 = 2wL j m{T) {d^f. The path integral of the 
Abelian gauge theory in the presence of a monopole with charge ia, located at position x is equivalent to 
the insertion of e ±lQj ' CT( - x ' into the path integral over the dual scalar fields [22, 27]. The complete partition 
function of the long distance effective theory is a sum over all topological sectors, each of which may contain 
an arbitrary number of monopoles and antimonopoles (whose charges sum to give the appropriate topolog- 
ical class). Summing over all numbers and locations of monopoles (and antimonopoles), weighted with the 
appropriate fugacity, directly yields the result (3.25). 
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interaction vertex has the additional shift symmetry 



cr -> cr + 2vr/ij, i=l,...,iV-l (3.26) 

where {//i} are the iV— 1 fundamental weights of the SU(N) algebra. These are defined by 
the reciprocity relation with the simple roots, 

Hi ■ ctj = \5ij a 2 = 5ij (3.27) 

for i = 1, • ■ • , N—l, which implies that the fundamental weights {/Xi} form a basis which is 
dual to the fundamental roots {<x,}. The presence of the symmetry (3.26) is related to the 
fact that the vacuum of the original theory can be probed by N—l different types of external 
charges, distinguished by their (non-zero) values of iV-afity. This will be discussed below. 
Including a non-zero theta parameter in the original Yang-Mills action, 

tYM . C YM C 



r^r + iO — -2 tr F^F^ , (3.28) 

has the effect, in the grand canonical partition function (3.22), of multiplying the monopole 
fugacity by e %e ^ N and antimonopole fugacity by e~ %e l N . In the dual representation (3.25), 
this amounts to shifting the argument of the cosine by 9/N, so that the interaction becomes 

N 

-(^2cos{ai-(T + i6/N). (3.29) 

i=l 

In this form, 2ir periodicity of the theory with respect to 9 is not manifest. However, a shift of 
the dual scalar fields, cr — > a + (9/N)(3 with (3 = (0, 1, 2, • • • , N—l), converts the interaction 
term to the manifestly 2tt periodic form 



C 



N-l 

cos(«j • cr) + cos(«Ar • cr + 9) 

i=l 



(3.30) 



in which theta dependence only appears in the term involving the affine root. For simplicity, 
in the following subsections we will focus on the case of 9 = 0. 

3.1 Mass gap 

The cosine potential in the dual action (3.25) generates a mass term for the photons. Rescaling 
cr to put the kinetic term into canonical form and expanding the potential to quadratic order 
around the minimum at cr = gives 

N 



V{oi) = (const.) + ^ (cr i+ i - dj) 2 , (3.31) 



with (Jn+i = fi and 



{27T? LC-Am* ( 2?fX 



m: * = ^J-L( = Ami v [^-\ e~ AS e -8- 2 /(^ 2 (-w-)) . ( 3 .32) 
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A Zjy Fourier transform, 

N 



LEe-^ffj, p = 0,---,N-l, (3.33) 

i=i 

diagonalizes this mass term and yields 

AT-l 

^(cJi) = (const.) + | m l \°p\ 2 ' ( 3 - 34 ) 
p=0 

with 

m p = m 7 sin ^— ■ . (3.35) 
Expressing m 7 in terms of the renormalization group invariant scale A, defined by 

A fe ° = /> {Ng 2 (ri)~ bl/b0 e^'/WO*)) (3. 36) 

(with &o = H/3 and b\ = 17/3), yields 

m-y = A A (AAL) 5 / 6 |ln AAL| 9/U , (3.37) 

where A is an 0(1) coefficient. Relative corrections suppressed by powers of g 2 {mw) ~ 
l/|lnAAL| have, of course, been neglected. The result (3.35), for p = l,--- ,N—1, shows 
that the A— 1 photons of the "unbroken" U(l) N ~ l gauge group receive non-zero masses due 
to nonperturbative effects. 15 

3.2 String tensions 

Let us first examine the vacuum structure of the dual theory in more detail. The dual scalars 
are defined to be periodic with period 2tt. This implies that shifting a by 2ir times any root 
vector is an identity, <x = a + 2irai for all a« G As noted earlier, the dual action (3.25) 
is also invariant under a — > cr + 2^^^, where {/^i} are the fundamental weights of the SU(N) 
gauge group, defined by the reciprocity relations (3.27). The simple roots {oti} generate the 
root lattice A r . Its dual, the weight lattice is generated by the fundamental weights 
The root lattice is a sublattice of the weight lattice and their quotient is 

A w /A r = Z N . (3.38) 

This implies that the dual theory potential, V(er) = — C cos(aj -<r), has A isolated minima 
lying within the unit cell of A r . These minima are located at cr = and 

a = 2nnj , j = 1, ■ ■ ■ , N-l . (3.39) 

(Equivalently, one may describe the minima as lying at 2tt j fi\, for j = 0,1,- •• , A, since 
= j A*i + a f° r some a G A r .) 



1 The vanishing p = mass corresponds to the extra decoupled photon which was added to simplify the 
duality transformation but is not present in the original theory. It should be ignored. 
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Let 1Z be some chosen irreducible representation of SU{N). The expectation value of 
the Wilson loop Wn(C) characterizes the response of the system to external test charges in 
the representation 1Z. In a confining phase with a non-zero mass gap, if external charges in 
representation 1Z cannot be screened by gluons, then expectation values of large Wilson loops 
in this representation are expected to decrease exponentially with the area of the minimal 
spanning surface, 

(W n (C))^e~ T ^ Arc ^ . (3.40) 

Here S denotes the minimal surface with boundary C, and T{7Z) is the string tension for 
representation 1Z. Such area law behavior implies the presence of an asymptotically linear 
confining potential between static charges in representation 1Z and anti-charges in represen- 
tation 1Z, Vft(x) ~ T(TZ) |x| as |x| — > oo. 

The irreducible representation 1Z may be associated with its highest weight vector w £ 
A^. Identifying weight vectors which differ by elements of the root lattice produces a % 
grading of representations which corresponds to their A^-ality (the charge of the representation 
under the Z_/v center). In particular, if this equivalence associates the highest weight vector 
w with k times the fundamental weight m, 

w = k fj>i + q , for some a £ A r , (3-41) 

then the representation 1Z has A^-ality A;. 16 

As discussed in Refs. [27, 22], the insertion of a Wilson loop Wn(C) in a representation 
1Z with non-zero AT-ality k corresponds, in the low-energy dual theory, to the requirement 
that the dual scalar fields have non-trivial monodromy, 

/ d(T = 2TTH k , (3.42) 
Jc 

where C is any closed curve whose linking number with C is one. In other words, in the 
presence of the Wilson loop W-ji(C) the dual scalar fields must have a discontinuity of 2irnk 
across some surface £ which spans the loop C. One way to see this is to go back to the 
duality relation. For simplicity, consider the case of a large planar loop lying in the xy- 
plane. As the size of the loop grows, the spanning surface £ approaches an infinite flat plane. 
In the presence of the Wilson loop, the Abelian duality relation F ~ *da is replaced by 
F ~ *da + Hk 5(z) dx A dy. Therefore the dual scalars a must be discontinuous across T, in 
order for the field strength F to be continuous. 

The fact that dual low energy theory depends on the representation 1Z of the Wilson loop 
only through its A r -ality k shows that there are only N—\ distinct string tensions, referred 
to as fc-string tensions, {T^}. (Charge conjugation symmetry implies that = T/v-fc.) The 
dual theory representation of Wilson loops also shows that external charges in representations 
with zero A^-ality will not be confined. These are precisely the representations which can be 
screened by adjoint representation gluons. 



Representations contained in the product of m powers of the fundamental representation with n powers 
of the antifundamental have iV-ality m—n. 
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To evaluate a Wilson loop expectation value, one must minimize the dual action in the 
space of field configurations satisfying the monodromy condition (3.42). To extract the string 
tension, 

area(S)^oo area(2j) 

it is sufficient to consider the limit where S fills the xy-plane. In this case, the field <x(x) 
will only depend on z. It must approach some minimum of the dual potential at infinity, 
lim z - t ± OQ (t(z) = 2irfj,i, and must be discontinuous across z = with a jump given by the 
prescribed fundamental weight, lim 2 ^ + ct{z) — lim 2 ^ - a { z ) = 2TVfik (mod 2ir). Because 
shifts by 27r/Zfc are an invariance of the dual potential V(<r), one may equally well minimize 
the action for field configurations cr(z) which are continuous but whose asymptotic values 
differ, 

AS{a) 



Tfc = min 



(3.44) 

A.cr=2-KHf 1 (mod 2ir) 



r(z) area(I 

where Aa = cr(oo)— <x(— oo), and AS (a) is the dual action minus its vacuum value. Explicitly, 

2 



mm 
tr{z) 



. (3.45) 

Acr=2iT(j,i f . (mod 2tt) 



In other words, the fc-string tension T\. equals the action of a kink solution with topological 
charge k in this one dimensional theory. 

The width of the kink solution must be of order of the inverse photon mass m~ l . Con- 
sequently, the /c-string tension will have the form = T, where 

T = C/m 7 ~ A 2 (ALJV)" 1 / 6 |log(ALiV)r 3/11 , (3.46) 

and fk is an 0(1) coefficient. Even without finding the minimizing kink solutions explicitly, 
it is apparent that the resulting /c-string tension Tf. will be non-zero (for k = 1, • • • , TV— 1), 
and must satisfy the convexity relation T^+i < + T\. 

We were unable to solve the kink equations of motions analytically for general N, but 
when N = 2 the equations of motion reduce to Sine-Gordon model. In this case, one finds 

T = T 1 =4v / 2CM 7 . (3.47) 

3.3 Larger size or larger N 

The above semiclassical analysis of the deformed Yang-Mills theory is reliable provided there 
is a parametrically large separation of scales between the lightest VF-boson mass, m\y = 
2n/(NL), and the nonperturbatively induced dual photon mass m 7 . Their ratio scales as 



mw 



(LiVA)- 11 / 6 | log(Li\TA)|- 9 / u , (3.48) 



and hence there is a large separation of mass scales provided LNA <C 1. In this regime, the 
monopole gas is highly dilute and a semiclassical analysis is justified. Increasing LNA, by 
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increasing N, L, or both, decreases the separation of scales; the heaviest photon mass, m 7 , 
grows while the lightest W mass, mw-i drops. When LNA ~ 1, the scale separation is entirely 
lost, the effective 't Hooft coupling A = g 2 N at the scale of mw ceases to be small, and the 
long distance dynamics can no longer be described by a weakly coupled {/(l)^ -1 effective 
theory. 

One can consider sending N to infinity while staying within the analytically tractable 
regime. This is a double scaling limit in which g 2 N and LNA are both are held fixed (and 
both are much less than unity) as N — > oo. Taking a large N limit in this fashion allows 
monopole effects to survive and to continue dictating the nonperturbative physics of the 
deformed Yang- Mills theory. However, this region shrinks to a vanishingly small window in 
the large N limit, since the double scaling implies that < L -C L max with L max A ~ 1/N. 
For any fixed compactification size L, if one sends N — > oo the deformed YM theory ceases 
to possess a monopole dominated, Abelian long distance regime. 1 

3.4 Connection to integrable Toda theory 

The 7*n symmetric model (3.25) is a deformation of a complex affine Toda theory with action 

N 



S 



affine Toda 



I ^(£)W-C|>-<— •«>]• (3.49) 



This complex (CPT-noninvariant) action describes a plasma which is composed solely of 
monopoles with no antimonopoles. (Due to the existence of the affine root, one can have 
a neutral plasma composed solely of monopoles!) Interestingly, the soliton spectrum of the 
affine Toda theory is exactly computable. When reduced to one dimension, this theory is an 
integrable system as shown by Hollowood [40], using techniques due to Hirota [41, 42]. 

As discussed above, the £;-string tension Tk is equal to the action of the kink solution 
with topological charge k. Borrowing the exact soliton spectrum from Ref. [40], one finds 
that the /c-string tensions in the affine Toda theory are given by 

T affine Toda = TN ^ k = 1, N—l , (3.50) 

with T given above in Eq. (3.46). 

The long distance effective theory (3.25) for our deformed Yang-Mills theory (when 
LNA <C 1) is a deformation of the affine Toda system by complex conjugation. Unfortu- 
nately, unlike the integrable affine Toda system, when N > 2 the resulting CPT invariant 

17 An analog of this double scaled limit was previously discussed by Douglas and Shenker in mass deformed 
SU (N) Af — 2 supersymmetric Yang-Mills theory on R 4 [26] . This deformed Af — 1 supersymmetric theory, 
just like our deformed Yang-Mills theory, possess a regime in which the long distance gauge dynamics reduces 
to the Abelian subgroup l^l)^ -1 . Ref. [26] shows that in the N -> oo limit of the mass deformed theory, 
the Abelian long distance regime is preserved only if the mass deformation m is sent to zero in a correlated 
fashion, m/A ~ 1/N 4 . In particular, at any fixed non-zero m, if one takes N — > oo first then there is no 
regime of the supersymmetric gauge theory in which the long distance dynamics remains Abelian. Although 
this phenomena only appeared previously in the context of supersymmetric gauge theories, it is generic in 
deformed Yang-Mills and other deformed QCD-like theories. 
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system is no longer exactly integrable according to Hirota's criteria. Consequently, we do 
not expect fc-string tensions in the deformed Yang-Mills theory to have the sine-law form 
(3.50). 

Recently, there have been attempts [43], to model the strongly coupled confined regime 
of Yang-Mills theory assuming the Wilson line has the center-symmetric form (3.1). (See also 
earlier related work in Refs. [44, 45, 46, 39].) A few remarks concerning the connection with 
Ref. [43] may be in order. First, our results for the fe-string tensions do not support the claim 
of Ref. [43], which asserts that /c-string tensions will have the sine law form (3.50). As just 
noted, sine-law string tensions are a property of the affine Toda subsystem, whereas the center- 
stabilized Yang-Mills in a weak coupling regime is dual to a real deformation of the affine 
Toda theory. We see no reason to believe that the /c-dependence of the string tensions will be 
unaffected by the deformation. Secondly, it should be emphasized that the deformation (2.5) 
stabilizes the center symmetric vacuum in the weakly coupled regime, and thereby provides 
a window in which a semiclassical analysis is reliable. Many earlier discussions of center 
symmetric backgrounds do not clearly distinguish the weakly coupled "Higgs" regime, in 
which fluctuations of the Wilson line eigenvalues are small, from the strong coupling regime 
in which the eigenvalues have large fluctuations and are essentially randomized over the unit 
circle. In our deformed Yang-Mills theory, both regimes exist. As the compactification size L 
increases, the theory moves from the weakly coupled regime to the strongly coupled regime. 
These two regimes are expected to be smoothly connected — no physical order parameter 
sharply distinguishes the two regimes. Nevertheless, the long distance physics of the weak 
coupling Higgs regime is effectively Abelian, while in the strong coupling regime there is no 
length scale beyond which the dynamics can be described accurately in terms of Abelian 
degrees of freedom. 
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In the absence of the complex conjugate term in the potential, there is a change of variables which converts 
the soliton equation of motion into "Hirota bilinear type," which is synonymous with solvability [40]. The 
presence of the complex conjugate term spoils the bi-linearity. 



A. Monopole measure 



The appropriate one-loop measure for integrating over configurations of a single monopole 
(of any type) may be expressed as 19 

^monopole = ^ e^ 3 ™^ J a J ^ [det'(-^adj)] , (A.l) 

where a £ M 3 is the monopole position, (j) £ [— vr, ir] is the internal U(l) angle of the monopole, 
and fi is the (Pauli-Villars) renormalization scale. Global U(l) gauge transformations (in the 
U(l) subgroup associated with the given type of monopole) shift the angle <j>. Fluctuations in 
the position and U{\) angle of the monopole represent the four zero modes in the monopole 
small fluctuation operator; 20 the factor of fi 4 can be viewed as the contributions of the Pauli- 
Villars regulator fields associated with these bosonic zero modes. The exponential factor 
is, of course, the exponential of minus the classical action of the monopole. The collective 
coordinate Jacobians are given by [39] 

2irS 1/2 

J a = S YM , Jq = — = NL S YM , (A. 2) 

a j ■ A 

where A = {0,27r,47r, • • • , 2ir(N— l)}/(iVX) are the eigenvalues of — iln$7. The primed deter- 
minant represents the result of Gaussian integrals over all fluctuations other than zero modes; 
the prime on the determinant denotes omission of the zero modes. The contributions from 
gauge bosons and ghosts, 

[det'(-D\ v - 2iV) adj ] ~ 1/2 x det(D 2 ) adj , (A.3) 

V v ' V V ' 

gauge bosons ghosts 

combine to give this simple form because 

[det'(-^V " 2iV) adj r 1/2 = [det(- J D 2 ) adj ]" 2 (A.4) 

in any self-dual background. These functional determinants may be regularized using the 
Pauli-Villars scheme. 

The fields of fundamental monopoles reside entirely within an SU(2) subgroup of SU(N), 
and the characteristic size of these monopoles is given by the inverse of the lightest W- 
boson mass, tu^} ~ NL. (This is the only scale which appears in the classical equations 



19 The following summary is an adaptation of the appendix of Ref. [39], which treats the monopole measure 
in supersymmetric Yang-Mills theory. 

20 For comparison, recall that an instanton in SU(N) Yang-Mills theory has iN zero-modes. (See, for 
example, Ref. [47].) For SU(2) gauge theory, these are the four zero modes corresponding to changes in the 
instanton position, one for its size, and three global gauge rotations. For SU(N) with N > 2, there are 
in addition 4A?"— 8 zero modes associated with changes in the embedding of SU(2) within SU(N). When 
compactified in one dimension, with non-trivial holonomy f2, one may regard an instanton as being composed 
of Af independent monopole constituents [45, 44], each of which carries four zero-modes. 
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for the monopole.) The regularized scalar determinant depends on the cube root of the 
renormalization scale, det(— D 2 ) ~ /i 1 ^ 3 . Since the determinant is dimensionless, it must 
have the form 

[det(-L) 2 )]- 1 = 2vrC (/iNL)-^ 3 , (A.5) 

where C is a pure number (iV-independent). Consequently, the one-loop monopole measure 
equals 

^monopole = C^ 3 (NL) I (5 YM ) 2 e~ S ™ +AS d 3 d dtt . (A.6) 

Performing the trivial integral over the angle O, the result is the £ d 3 a, with £ the monopole 
fugacity. Choosing to use mw as the value of the renormalization point yields the expression 
(3.23) for the fugacity. 
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